Mon. Not. R. Astron. Soc. 000. HHT51 ('20031 Printed 2 February 2008 (MN IMpX style file v2.2) 



Massive black hole seeds from low angular momentum 
material 



Savvas M. Koushiappas, 1 James S. Bullock, 2 * Avishai Dekel, 3 

1 Department of Physics, The Ohio State University, 174 W. 18th Ave, Columbus, OH 43210 USA; smkoush@mps.ohio-state.edu 

2 Harvard Smithsonian Center for Astrophysics, 60 Garden St. Cambridge MA 02138 USA; jbullock@cfa.harvard.edu 

3 Racah Institute of Physics, The Hebrew University, Jerusalem, Israel; dekel@phys.huji.ac.il 



Released 2002 Xxxxx XX 



ABSTRACT 

We present a cosmologically-motivated model in which the seeds of supermassive 
black holes form out of the lowest angular momentum gas in proto-galaxies at high 
redshift. We show that under reasonable assumptions, this leads naturally to a cor- 
relation between black hole masses and spheroid properties, as observed today. We 
assume that the gas in early-forming, rare-peak haloes has a distribution of specific 
angular momentum similar to that derived for the dark matter in cosmological N-body 
simulations. This distribution has a significant low angular momentum tail, which im- 
plies that every proto-galaxy should contain gas that ends up in a high-density disc. In 
haloes more massive than a critical threshold of ~ 7 x 10 7 Mq at z ~ 15, the discs are 
gravitationally unstable, and experience an efficient Lin-Pringle viscosity that trans- 
fers angular momentum outward and allows mass inflow. We assume that this process 
continues until the first massive stars disrupt the disc. The seed black holes created 
in this manner have a characteristic mass of ~ 10 5 M Q , roughly independent of the 
redshift of formation. This serves as a lower bound for black-hole masses at galactic 
centers today. The comoving density of mass in black hole seeds grows with time, 
tracking the continuous production of critical mass haloes, and saturates when cos- 
mic reionization acts to prevent gas cooling in these low-mass systems. By z ~ 15, 
the comoving mass density becomes comparable to that inferred from observations, 
with room for appropriate additional luminous growth during a later quasar accretion 
phase. The hierarchical merger process after z ~ 15 naturally leads to a linear corre- 
lation between black-hole mass and stellar spheroid mass, with negligible black hole 
masses in disc-dominated galaxies. The formation of massive seeds at high redshifts, 
and the relatively important role of mergers in the buildup of today's black holes, are 
key elements in the proposed scenario. 
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1 INTRODUCTION 

Detailed studies of gas and stellar kinematics near the 
centers of present-day galaxies have revealed that al- 
most all galaxy spheroids (ellipticals and bulges of spi- 
rals) host super-massive black holes, with masses M^h = 
(10 6 - 10 9 )M o , or ~ 10~ 4 of the total stellar mass of 

their parent galaxi es fe. g. iKo rmendy &, Gebh ardd 120011 

iFerrarese fc MerritdliuOol: iMaeorrian et al.lll99gft "Although 
most of these black holes are not associated with quasar 
activity today, presumably they represent the now-dormant 
counterparts to the quasar-powering engines known to exist 
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at high redshift llFan et alJl200llh Indeed, under reasonable 
assumptions, the luminosity function of Active Galactic Nu- 
clei (AGN) can be explained by modeling mass accretion 
onto black holes that have the sam e range of masses as ob- 
served in local spheroids (see e.g. ISteed fc Weinbe'rel |2003| 
and references therein). The standard expectation is that 
once a primordial galaxy is populated with a "seed" black 
hole at some early time, the black hole can grow via the 
accretion of available gas and, in the process, give off light 
in proportion to its Eddington luminosity. 

The seed population is an important ingredient in any 
model that aims to explain AGN or local black holes, but 
the origin and nature of these objects remain unknown. In 
this paper, we provide a model for producing rather mas- 
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sive (~ 1O 4_5 M0 ) seed black holes at early times from high 
density, low angular momentum gas in rare halos. We argue 
that massive, early forming seeds help explain the high red- 
shift AGN population, and, unlike lower mass seeds, may 
give rise to the observed correlations between host galaxy 
properties and black hole masses without the need for fine- 
tuned feedback. 

Consider, for example the analysis of iBarth et alJ 

(2003), who have estimated a mass of Mbh = (2 — 6) x 
1O 9 M0 for a black hole powering a z = 6.4 quasar found 
in the Sloan Digital Sky Survey. For the ACDM cosmology 
adopted below, the age of the universe at this redshift is 840 
Myr. Even the remnant of a very massive, early-forming star 
of ~ IOOM0 would have to accrete mass at the Eddington 
limit for 17 e-folding times in order to grow to the observed 
mass. For a standard, thin-disk accretion model with an ef- 
ficiency ~ 0.1, the e-folding Salpeter time is ~ 45 Myrs, so 
an intermediate-mass seed would require 770 Myrs of ac- 
cretion 1 in order to reach to observed mass, corresponding 
to a formation redshift for the seed of Zf£, 40. This is an 
extremely early redshift in the context of ACDM structure 
formation, since the minimum halo mass for H2 cooling at 
this time represents a > 5a fluctuation. While this does not 
completely rule out the possibility that the seed was small, 
it does require an extremely rare event. A 1O 5 M0 seed, on 
the other hand, would require ~ 10 e-folding times of accre- 
tion, and thus a more recent formation redshift zf^i 11. Of 
course these numbers are sensitive to certain assumptions, 
and may be influenced by allowing low efficiency accretion 
or super-Eddington luminosities. Nevertheless, it is gener- 
ally easier to understand the population of z> 6 quasars if 
massive seeds are present (see also Haiman & Loeb 2001 for 
a related discussion). 

Perhaps the most intriguing clues to the nature of black 
hole formation and growth are the strong correlations be- 
tween black hole masses and global properties of their host 
spheroids — scales that differ by many orders of magnitude. 
The black hole mass is tightly correlated with the central 
velocity dispersion of the spheroid, and shows a significant 
correlation with the spheroid luminosity or mass. The theo- 
retical challenge is to connect the galaxy-formation process 
originating on scales of megaparsecs with the black-hole for- 
mation processes on sub-parsec scales. In addition, the epoch 
of spheroid formation is likely much later than the time when 
the first quasar black holes begin to shine, so an additional 
challenge is to connect these two components in time, in 
a way that maintains, or even creates, the observed local 
correlations. 

Several scenarios have been presented for the for- 
mation of massive black hole seeds. The general frame- 
works range from those that rely on the coalescence 
of massi ve stellar remnants tha t sink to the galaxy 
cent ers jpuinlan fc Shapiro! 1199(1 Iffills fc Bendeil Il993 
iLed Il995l : lEbisuzaki et al.1 1200 J) to one that models 
black hole growth via ballistic p articl e capt ure dur- 
ing bulge collapse jAdams et al.1 l200ll 120031) . Other 
models are based on the hydro-dynamical evolution of 
super-massive objects formed directly out of primor- 



1 This timescale is much lar ger than the e stimated lifetimes of 
quasars at lower redshift (see lMartinill20o3 for a recent review). 



dial gas jHaehnelt fc Reed Il993t lLoeb fc Rasid Il994l; 
[ Eisenstein fc Loebl Il995al: Ha ehnelt. NatamHan fc Reed 
Il998t lGnedinll200lt iBromm fc Loebl I2003D . Our model is 
related most closely to this third type of scenario. 

The key obstacle in any model of black hole formation 
is clearly the centrifugal barrier. In the standard cosmolog- 
ical framework of hierarchical structure formation, proto- 
haloes acquire angular momentum by tidal torque s exerted 
by the background fluctuation field llHovldll949t I Feeble: 



19691: iBarnes fc EfstathiorJll987tlPorciani. Delel fc HoffmaH 



2002iS~ Haloes obtain a specific angular momentum that 
is on average 10 — 100 times smaller than that required for 
centrifugal support at the virial halo size (with a spin pa- 
rameter of A ~ 0.01 — 0.1, see below). This is roughly inde- 
pendent of mass and time. Thus, even after cooling, most of 
the gas can contract only to radii 10-100 times smaller than 
the initial halo radius before it is halted in circular orbits. 
The implied scale (at z — 0) is ~kpc, many orders of mag- 
nitude larger than the corresponding Schwarzschild radius, 
~ 1O _5 (M/1O 9 M0 )pc. The situation may get even worse 
at earlier times because the typical collapsing mass M*(a) 
is smaller in proportion to a high power of the expansion 
factor a = (1 + z)^ 1 . This implies that the ratio of the cen- 
trifugally supported size (or. aMi(fl) 1 ' 3 ) to the Schwarzschild 
radius (oc M*(a) if the black hole mass is a constant frac- 
tion of the halo mass) increases dramatically at early times, 
unless the sites of black hole formation are rare haloes much 
more massive than M*(a). This line of argument indicates 
that the seed black holes originate from material at the very 
lowest end of the angular-momentum distribution within 
rare, massive haloes at high redshift. This material should 
then perhaps lose all its angular momentum via an effective 
mechanism of angular momentum transfer, which would be 
enhanced i n such a high-densit y environment. 

lEisenstein fc Loebl (u995af) considered the sites of black 
hole formation to be the very high density discs which form 
in those rare haloes that have extremely low total angular 
momentum. Assuming a standard CDM scenario, they used 
an analytic model to predict the distribution of halo spins 
and estimated the abu ndance of haloes with low e nough spin 
for efficient viscosity feisenstein fc Loeblll99 5b). They ar- 
gued that there may be enough low-spin haloes to explain 
the black hole abundance. They then adopted the common 
view that viscosity in a differentially rotating gaseous disc is 
the mechanism that transfers angular momentum outward, 
thus allowing gas to contract into a black hole. In their cal- 
culations, they considered hydrodynamic (or perhaps hydro- 
magnetic) turbulence to be the viscous mechanism respon- 
sible for the angula r momentum loss, and ap pealed to the 
a-disc prescription JShakura fc Sunvaevlll973l) . 

Our approach is related, but qualitatively different. We 
propose a scenario where seed black holes form from mate- 
rial with low angular momentum in all haloes that are mas- 
sive enough to host an unstable self-gravitating disc. The 
required high-density discs can only occur in haloes that 
collapse at high redshift, z > 10 — 20, when the angular mo- 
mentum distribution of the (unprocessed) gas should be the 
most similar to that of the dark matter, and b efore reioniza- 
tion or some other feedb ack mechanism (e.g.. iDekel fc Silkl 
ll986tlDekel fc WodEioll) prevents the growth of cold discs. 
The crucial point is that every halo, regardless of its to- 
tal angular momentum, is expected to have a distribution 
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of specific angular momentum including a tail of low-spin 
material. 

The haloes we consider are massive enough so that the 
baryons trapped in their potential wells can cool and con- 
dense. In the absence of any other effect that may alter the 
angular momentum content of baryons, the lowest spin bary- 
onic material in each system ends in a high-density central 
region forming a proto-galactic disc. As an illustrative fidu- 
cial model, we assume that the angular-momentum distri- 
bution in each halo m atches the "universal" distribution of 
iBullock etall <2001al) (B01 hereafter). If the angular mo- 
mentum is approximately preserved during baryonic infall 
into the disk, the B01 angular momentum distribution im- 
plies that the inner disk becomes self gravitating with the 
surface density inversely proportional to the radius. We show 
below that if this behavior extends below the inner ~ 2% of 
the disc radius, then enough mass may acquire a high enough 
density for viscosity to efficiently redistribute the angular 
momentum of the disc. For the final contraction of disc mate- 
rial into a black hole we evaluate the angular-momentum loss 
by the effective kine matic viscosity that arises from gravi- 
tational instabilities l|Lin fc Pringlelll987li . This mechanism 
is somewhat more efficient than that from a typical a-disc 
model. This viscosity-driven mass inflow cannot continue in- 
definitely; once the disc significantly fragments, or once it is 
heated or even disrupted by feedback, the inflow stops. 

The semi-quantitative scenario presented here can be 
regarder as an illustration of a general model based on the 
assumption that proto-galactic gaseous discs develop high 
density inner regions dominated by low angular momentum 
material. Our current analysis is basically a feasibility study 
of such a model, and an attempt to recover its robust predic- 
tions and to constrain the model parameters by theoretical 
considerations as well as by the observations. 

The outline of this paper is as follows. In §|2|we discuss 
the angular momentum distribution in proto-galaxies and 
the implied surface density profile for the cold discs. In § |3 
we investigate the effect of kinetic viscosity on the evolution 
of the disc. In §|l]we calculate the resulting seed black hole 
mass and mass function. In § |S] we calculate the redshift 
evolution of the seed black hole mass density. In § |S|we make 
a preliminary attempt to address the correlation between 
the masses of black holes and galactic stellar spheroids, to 
be pursued in a subsequent paper. In § [7| we summarize 
our conclusions and discuss our results. When necessary, we 
assume a ACDM cosmology, with a Hubble constant h = 0.7 
(Ho = lOO/i km s -1 Mpc^" 1 ), universal mass density Qm = 
0.3, cosmological constant Qa = 0.7, and rms fluctuation 
amplitude ag = 0.9. 



2 SPIN DISTRIBUTION AND DISC PROFILE 

We focus on the proto-galaxies forming in the redshift range 
z ~ 10 — 30. In the standard picture of structure forma- 
tion, galaxies form as a uniform mixture of dark matter and 
baryons in small over-densities that grow by gravitational 
instability. Until its expansion turns around into a collapse, 
each halo acquires angular momentum through tidal torques 
induced by the density fluctuations around it. The dissipa- 
tionlcss collapse brings the dark matter into virial equilib- 
rium in an extended spheroidal halo. In the presence of an 



efficient cooling mechanism, gas cools and condenses to a 
centrifugally supported disc. 

Atomic line cooling can bring the gas to ~ 10 4 K, and 
further cooling down to ~ 300K requires the presence of 
molecular hydrogen. The minimum halo mass for atomic 
cooling is M a ~ 6 x 1O 7 M [(1 + z)/18]" 3/2 , and for molec- 
ular hydrogen survival it is Mh 2 ~ 5 x 10 6 M W [(l + z)/18l~ 7 
in the redshift range z ~ 13-30 llTegmark et all Il997l: 
lAbel et al]ll998T) . The haloes we consider below are just 
above the limit for atomic line cooling and significantly 
above the limit for molecular hydrogen cooling. We therefore 
expect the gas to cool and settle in an angular momentum 
supported disc on a timescale comparable to the dynamical 
time: 

t dyn = (G Pvir )- 1/2 ^77MyrA 178 (z)- 1/2 (I±£) , (1) 

where we have assumed our standard cosmological param- 
eters and Ai78(z) is the virial overdensity in units of 178 
times the background density. 2 . Hereafter, due to the small 
redshift dependence of A v i r at the redshifts of interest we 
set A 17 $(z) = 1. 

Consider a halo of virial mass M v and spin parameter 3 
A = J[\/~2M V V V R V \~ , where J is the total halo angular mo- 
mentum. The virial radius isi? v = 384 pcM^ 3 [(l+2)/18] _1 
and the virial velocity is V v = 11 kms -1 M 7 1/3 [(l + z)/18] 1/2 . 
Here, M 7 = M V /10 7 M Q . B01 found that in each halo the 
cumulative mass with specific angular momentum less than 
j is described well by the function 

M(<j)=fj,M v ^—, fj,>l. (2) 
3o+3 

This profile has an implicit maximum specific angular mo- 
mentum j max = jo/(/i— 1), where jo = V2V V R v X/ b(fi) , with 
fe(/x) = — /xln(l — /i -1 ) — 1. The parameter n characterizes the 
shape of the angular momentum distribution. When fi 1 
the deviation from a power law is pronounced with a larger 
fraction of the mass having low specific angular momentum, 
while if n 3> 1 the distribution resembles a pure power law. 
The key property of this distribution is that at least half 
of the halo mass is in the power-law regime, M(< j) oc j, 
and the simulations confirm that the power-law behavior ex- 
tends at least over two decades in mass, down to ~ 1% of 
the total halo mass. 

The B01 pr o file w as confirmed bv lChen fc Jind (12002ft 
and IColm et alJ ll2 002) using N-body s imulatio ns, and by 
Ivan den Bosch et alJ fcOOft.lchen. Jing fc Yoshikawl (1200311 
and ISharma fc Steinmetzl J2003T) using adiabatic hydrody- 
namic simulations. These studies reveal that the B01 angular 
momentum distribution remains a good fit to haloes at red- 
shifts as high as z = 3. In what follows we will assume that 

2 If we use the approximati on of spherical top-hat collapse given 
in iBrvan fc Normanl <1998l) we get A vir (z) ~ (18tt 2 + 82x - 
39x 2 )/n M (z), where x + 1 = n M (z) = f2 M (l + z) :i /[n M (l + 
z) s + S7a] and S2\i (2) is the ratio of the mean matter density to 
critical density at redshift z. In the ACDM cosmology that we 
adopt, A v i r (z = 0) ~ 337 and A v i r (z) — > 178 at high redshift, 
approaching the standard CDM (Qm = 1) value. 

3 This spin parameter has been proposed by B01 as a practical 
modification of the conventional spin parameter, defined as A = 
J^/|-B|/GMv //2 , where E is the halo energy. The two definitions 
of A coincide for a singular isothermal sphere truncated at R v . 
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Figure 1. The mass per unit specific angular momentum for 
different values of the parameter £. Note that in our fiducial case 
of £ = 1 (B01), this function approaches a constant at small j. 



the distribution holds at much higher redshifts, z ~ 15. Since 
the ACDM universe is effectively Einstein-de Sitter beyond 
z ~ 1, and the power spectrum is effectively a power law in 
the relevant mass range, one expects self-similarity between 
the haloes at z ~ 3 and at z ~ 15. Therefore, the profile at 
z ~ 15 can be assumed to resemble the B01 sha pe derived at 
z = 0. Simple theoretic al considerations (B01, iDekel et alJ 
l2000l:lMaller" et al . 2002) argue that this angular momentum 
distribution is a natural outcome of the hierarchical cluster- 
ing process. A qualitatively similar distribution is obtained 
when linear tidal torque theory is applied to haloes shell 
by shell, or when the orbital angular momentum of merging 
haloes turns into spin of the product halo via tidal stripping 
and dynamical fr iction. Other justifications follow from more 
detai led models jMaller fc Deke]|l2002t Ivan den Bosch et alJ 
l200ll) . These analytic models provide addivtional support to 
the assumption that a similar average halo profile is valid at 
low and high redshifts. 

As mentioned above, the B01 distribution is a good 
fit to simulation results for angular momentum as small as 
~ 1% the characteristic angular momentum of the halo (cor- 
responding to ~ 1% the characteristic disc radius. As dis- 
cussed below (see Eq. 13), our calculation assumes that the 
B01 form holds on scales below ~ 2% of the characteristic 
angular momentum, so only a mild extrapolation of the low 
angular momentum tail is required, and being a power law, 
this extrapolation seems plausible. Nevertheless, the angular 
momentum distribution of the gas may not exactly mirror 
that of the dark matter. In order to evaluate the qualita- 
tive effects of this possibility and the explore alternative 
extrapolations, we will also condier deviations from the B01 
distribution. 

U nlike the analysis in B01 ( but see |Bjillock_et_alJ 
| 2002l). other work s suc h as IChen fc Jind feOOSlC 
lyanden Bosch et alJ i2002T) and IChen. Jing fc Yosbdkawl 
l2003Ti ~focused on estimating the fraction of dark mass 
with negative angular momentum (in the direction of the 



total halo spin), and claimed that this fraction may be 
as large as ~ 20%. 4 By using the B01 profile as is we 
actually adopt the conservative approach which minimizes 
the efficency of black-hole formation, as follows. If a fraction 
of the gas starts with negative spin, it will mix dissipatively 
with the positive-spin gas leading to a purely rotating 
disc with no negative component. The implications on the 
disc formation depends on the way by which this mixing 
occurs. In the limiting case where the mixing occurs after 
the gas has already condensed to the disc, the negative 
component will mix positive spin m aterial of exactly the 
same distribution of spin magnitudes, llvan den Bosch et alJ 

2002) . This results in a fraction of mass with near-zero 
angular momentum, which would increase the efficiency of 
black-hole production. If angular momentum mixing occurs 
during the process of gas cooling and collapse before the 
disc forms, there is no a priori reason for the negative 
angular momentum material to mix with only positive 
angular momentum material of the same magnitude. If this 
mixing is uniform across the angular momentum profile, 
then the resulting distribution should be similar to the B01 
form, but with a slightly smaller overall spin parameter. 
The B01 form, ignoring the negative tail, is therefore the 
conservative case of minimum blak-hole formation. 

We assume that the specific angular-momentum distri- 
bution of the baryons resembles that of the dark matter 
halo. This assumption has been found to hold reasonably 
well in b oth radiative and non-radiative hydrodynamic simu- 
lation s, llvan den Bosch et alJl2002l IChen. Jing fc Yoshikawl 

2003) . This assumption should be valid especially in sys- 
tems that collapse early, as baryons would not have been 
affected by feedback associated with star formation. In ad- 
dition we assume that the angular momentum of gas is con- 
served as gas cools and falls to the center forming a centrifu- 
gally supported disc. This assumption is backed by obser- 
vations, where disc sizes seem to indicate that the angular 
momentum of the gas gained by tidal torques is retained. In 
addition, B01 found that the direction of angular momen- 
tum is well aligned at different radii, suggesting that the 
formation of a disc is indeed possible. Therefore, if there are 
no significant changes in the specific angular momentum of 
the baryons during the formation of the disc, then the disc 
angular-momentum distribution is Md(< j) = foM(< j), 
where fo is the final mass fraction in cold baryons. For the 
circular orbits in the disc one has j(r) = rV(r), so one can 
map the Ma(< j) distribution to a radial mass profile in the 
disc. Note that if Md(< j) oc j (as with the B01 distribu- 
tion at small j) and the disc has a self-gravitating rotation 
curve V 2 (r) ~ GMa(r)/r, then V(r) =const. and Md(r) oc r 
regardless of the initial properties of the halo and the pos- 
sible effects of baryonic contraction on it. This implies a 
disc surface density profile of S(r) oc r _1 in the inner disc. 
Since we expect these general trends to hold regardless of 
our assumptions about the halo profile, we make the simple 
approximation that the final total mass profile of halo plus 
disc follows that of an isothermal sphere, with M(r) oc r 



4 See the extensive discussion in B01 on the pitfalls that may arise 
when attempting to measure the negative fraction of material in 
dark halos. 
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Figure 2. The disc surface density profile [eq. JHJ] for different 
values of £. Our fiducial case, £ = 1, shows a steep increase toward 
small radii, while the £ = 2 case has a flat core, similar to an 
exponential disc. 



and V(r) = V v . In this case, the disc mass profile implied 
by the distribution is 



M d (f) = fnM v -!—, r = r/r d , 
1 + r 

where r d is a characteristic disc radius, 



(3) 



V2XR V 
6(/i) 



22pc Ao.. 



^(Mrx 173 ^)" 1 (4) 



and Ao.o4 = A/0.04. At r <C »*d, we have again Md oc r. This 
profile is valid out to a maximum radius f max = l/(p — 1)- 
The constant / is the fraction of cold baryonic mass in the 
halo that has fallen into the disc by redshift z (/ ^ /o). 

Recall that it takes a finite time, of order ~ tdyn, for 
a centrifugally supported disc to completely form after the 
baryons cool. Thus / in eq. JHJ becomes equal to the total 
cold mass fraction /o only after a time t ^ td yn - Before that 
time, we assume / = /o (t/tdyn), such that the disc mass 
grows linearly with time until t — td yn , 



M d (t) = fM v 



f M v [ — 

t'dyn 



t tdy 



(5) 



Hereafter we assume that the total cold mass fraction is 
fo = 0.03, which is roughly ~ 1/5 of the universal baryon 
fraction of the s tandard ACDM Un iverse we assume in this 
paper (see also iKlypin et ail 120021) . This is a conservative 
estimate for the fract ion of baryons that can form enough H 2 
to cool significantly dTeemark et aJll997tlAbel et al.lll99ct . 

Since our approach relies on an uncertain extrapola- 
tion of the B01 distribution, we explore a more general form 
for the angular momentum distribution, of which B01 is a 
special case. We consider the family of angular-momentum 
distributions that correspond to the following disc surface 
density profiles, parameterized by f , 



Sd(f) 



NffiM v £f 



«-2 



2-ktI 



(l + f«) 2 



r<(/t-l) : 



(6) 



with the normalization constant N = — l) 4 + 1]. 

In Fig. 1 we show the differential angular momentum 
distributions for different values of £. The specific case of 
f = 1 is the B01 form, while larger (smaller) values of £ 
correspond to distributions with less (more) low angular 
momentum material. In Fig. 2 we show the resulting disc 
surface density profiles for different values of £. Note the 
steep increase in density for f = 1 (B01), while for £ = 2 the 
density profiles are similar to an exponential disc. 

Throughout this paper, we focus on the regime where 
r< 1, and the surface density is well-approximated by 



Sd(r) 



(7) 



We assume the B01 limit for the most part of this paper 
(which corresponds to £ = 1), but also investigate how our 
results change with varying the value of £. 



3 ANGULAR MOMENTUM LOSS 

After baryons cool and settle into a dense inner disc, an- 
gular momentum is redistributed as a result of viscosity. 
We adopt the kinem atic v iscosity prescription propos ed by 
iLin fc Pringlel lll987l) (see iLauehlin fc Rozvczkalll996l for a 
detailed test of this picture). According to this mechanism, 
global gravitational instabilities in the disc give rise to a 
transfer of angular momentum outward, which can be mod- 
eled as a local kinematic viscosity. 

Consider a disc with local surface density E, angular 
frequency u>, and temperature T, expressed as a sound spe ed 
c a ~ 2(T/300K) 1/2 kms~ 1 . As discussed bv lToomreHl964h . a 
disc is gravitationally unstable as long as the minimum size 
of shear-stabilizing disturbances L s ~ G Elj -2 , is larger than 
the maximum size of dispersion-supported regions, Ld — 
CsG _1 E _1 . This condition is expressed by evaluating the 



Toomre parameter Q, 



Q 



Ld 



1/2 



GE' 



(8) 



In discs with Q < 1, the interaction of the global instability 
perturbations with the background shear field drives a mass 
inflow. The transfer of angular momentum occurs across re- 
gions of maximum size ~ L s over a typical timescale ^ uj^ 1 
due to an effective viscosity 

(9) 



Li 



G E 



In what follows, we use ui 2 — GM(r)r~ 3 . 

A typical viscous time for mass inflow at a radius r in 
the disc is t„ ~ r 2 v~ x . As mass flows in, the central re- 
gions of the disc become denser, so the value of Q decreases. 
If at any point Q became as small as Q 2 < H/r, where 
H ~ c s /u) is the height of the disc, then the instabilities 
would have become fully dynamical, fragmenting the disc 
into self-gravitating bodies, which would have stopped the 
angular momentum transfer. However, this effect is likely to 
be unimportant over the timescales of concern here because 
while Q is decr easing, the ratio H/r is actu ally decreasing 
at a faster rate iMineshige fc Umemurall997l) . thus delaying 
the global fragmentation and allowing the viscosity to act 
until some external process disrupts the disc. The underly- 
ing assumption is that while the global instabilities give rise 
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to density waves, which cause the transfer of angular mo- 
mentum, the disc does not fragment instantaneously. Star 
formation can occur concurrently with the viscous process 
until feedback from the first massive stars stops the gas in- 
flow either by heating the disc or by disrupting it altogether 
(if its virial velocity is small enough). 



3.1 Critical halo mass for viscous transport 

Only haloes more massive than a critical minimum mass 
host unstable discs in which the angular-momentum loss is 
efficient and a seed black hole can form. This is because, 
from eq. JHJ, we expect Q ~ M~ 1//3 , implying a critical 
mass above which Q < 1. For the specific case of the B01 



profile the scaling is Q oc M c 



-1/2 1/2 ,,- 

r d < oc M d 



1/3 



(1 + *) 



-1/2 



independent of radius. If we assume Md $S foM v , then set- 
ting £ = 1 in eq. yields a lower limit for host halo masses 
hosting Q < 1 discs: 



M':. 



> 1.2 x 1O 7 M T^ 2 



3/2 



A 



0.046 (n) 



0.03 

fo 



1.25 



3/2 



(10) 



All the quantities in brackets are of order unity, and T300 = 
T/300K is the typical temperature for gas that has cooled 
via molecular Hydrogen. 

A slightly improved estimate of this critical halo mass 
can be obtained by considering potential disc perturbations 
by major mergers during the gradual growth of disc mass 
until it reaches the baryon fraction fo at t ~ tdyn after viri- 
alization. If we approximate the rate of mass growth in the 
disc by eq. ©, then Q oc M^ 1 ^ 3 oc t _1//3 implies that ini- 
tially self-gravity is negligible compared to gas pressure. As 
more material is accreted, Q decreases, and it may eventu- 
ally reach Q 1 where angular momentum transfer begins. 
We define the timescale £q as the timescale needed for the 
disc to become gravitationally unstable had the accretion 
been uninterrupted. To evaluate this time, we assume that 
the disc is growing according to eq. @, compute the surface 
density from the B01 profile (eq. Q with £ = 1), and sub- 
stitute in eq. (JHJ for Q. We then set Q = 1 for the stability 
threshold and obtain 



t Q = 1.4tdy»T 3 ooM. 



-2/3 



(ifi) 



0.03 

fo 



1.25 
i 1 



A 



0.04%) 



(11) 



Figure 3 shows tQ as a function of host halo mass at 
z = 17 and at z = 13. We assumed the fiducial values 
£ = 1, T = 300.ff, f = 0.03, n = 1.25, and A = 0.04. 
These curves are compared to the virial dynamical times 
tdyn at the two corresponding redshifts. The intersection of 
the tQ and td yn curve marks the uninterrupted critical mass 
at that redshift, as estimated in eq. IjlOfl . However, a ma- 
jor merger may destroy the disc and disrupt the gas in-fall. 
If the typical time for a major merger to occur after the 
halo virializes is t mm , the disc reaches a maximum mass of 
Md (tmm/idyn)/oM v < foM v . In order to obtain a given 
critical disc mass where viscosity efficiently redistributes an- 
gular momentum, we require £q < t mm , i.e., the disc should 
become self-gravitating before it gets disrupted. In this case, 
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Figure 3. Timescales as a function of halo mass. Shown are (a) 
the time needed for the disc to become self-gravitating (tQ), (b) 
the time between virialization and the following major merger 
(tmm), and (c) the dynamical time (tdyn). Thick (red) lines cor- 
respond to z = 17 and thin (blue) lines to z = 13. The critical 
mass for a halo that can form a seed black hole is determined by 
the intersection of the corresponding lines for tQ and t mm . 



eq. il l 01 is an underestimate of the critical halo mass needed 
for disc instability that leads to efficient angular momentum 
transfer. In the regime we consider, tmm is always less than 

^dyn ■ 

Figure 3 also shows the characteristic timescale for ma- 
jor mergers t mm , in the extended Press-Schechter formalism 
JPress fc SchechteJ [1973 iBond et all Il99li lLacev fc Colel 
1993). Specifically, we associate tmm with the peak of the 
probability distribution for the time it would take a halo t o 
double its mass (see equation 2.21 in lLacev fc Cole!ll993r) . 
The intersection of the tQ and t mm curves (for a given red- 
shift) marks the more realistic critical minimum mass, which 
turns out to be roughly a factor of ~ 4 larger than the un- 
interrupted estimate (eg. HOI . 

Figuer 4 shows M£ calculated in the way just de- 
scribed as a function of redshift. The thick solid line was 
computed assuming the B01 profile for the gas density dis- 
tribution. The redshift scaling is similar to that estimated 
in eq. II10L confirming that the critical mass is varying only 
slowly with redshift. The critical mass is about a factor of 
~ 4 larger than the uninterrupted estimate throughout the 
whole redshift range of interest. Shown for reference is the 
minimum halo mass for molecular hydrogen to survive as 
an efficient coolant; it justifies our assumption that gas can 
cool in critical-mass haloes. 

Also shown in Fig. 4 is our estimated critical mass 
for a flatter surface density profile (or narrower angular- 
momentum distribution) with £ = 1.3. As would be ex- 
pected, M" A is higher here than it is for the £ = 1 case. 
This is because densities at small radii are lower for a fixed 
mass disc, so a more massive host is needed in order to 
supply the extra disc mass required for Q < 1. Unfortu- 
nately, the instability calculation is less straightforward for 
£ ^ 1 because Q is not the same at all radii in the disc: 
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Figure 4. The minimum mass of a halo capable of forming a seed 
black hole as a function of redshift. The thick solid curve marked 
B01 refers to our fiducial £ = 1 distribution. The thin solid curve 
refers to £ = 1.3, with a longer time of 20 Myr assumed for star 
formation. The dotted lines correspond to 2-8 cr fluctuations in the 
random fluctuation field, illustrating that the minimum masses 
of concern are rare systems. The short dashed line represents the 
minimum mass for atomic line cooling only, while the long dashed 
line is an approximation for molecular hydrogen cooling . 



Q oc f f 1- ^/ 2 . Indeed, for £ > 1 discs become unstable from 
the outside in; however, once inflow begins, mass transport 
always tends to increase the central density with time and 
therefore it is natural to assume that the central disc will 
become unstable within a viscous timescale. We can crudely 
estimate the critical mass in this case by evaluating Q and 
the viscous timescale at some characteristic radius. If we as- 
sociate the viscous timescale with the timescale over which 
the disc would get disrupted (see next section), we can de- 
fine the critical halo mass for black hole seed formation to 
be the mass scale at which Q < 1 for the given viscous 
timescale. The £ = 1.3 curve in Fig. 4 has been computed 
assuming t = 20 Myr. Despite the fact that this result re- 
lies on very crude approximations, it does serve to illustrate 
how narrower angular momentum distributions require more 
massive (and less abundant) haloes as the sites of black-hole 
formation. We come back to the implications of this later in 
the paper. 



3.2 Viscous inflow 

Consider now haloes more massive than the critical mass 
for disc instability and assume that they have formed discs 
with surface density profiles described by eq. Q. The vis- 
cous process transports mass inward from a radius r over 
the viscous time t v ~ r 2 f _1 . We make the assumption that 
if viscosity has been acting for time t, then all of the disc 
material originally within r v i s = (ut) 1 ^ 2 has lost its angu- 
lar momentum. Assuming a viscosity given by eq. this 



radius can be written as 




where r v i s = r v i s /r d and / = /o(tQ/tdyn)- For the special 
case of £ = 1, we have 




We see that the inflow radius grows linearly with the time 
available for viscosity to efficiently transfer angular momen- 
tum from the inner parts of the disc to outer radii. 

According to the above equations, if the viscous 
angular-momentum transfer were to continue uninterrupted 
for a long enough time, then the whole disc would have even- 
tually lost its angular momentum and collapsed. However as 
mentioned before, the major merger timescale sets an upper 
limit on the time available for this process, and in fact the 
feedback from star formation is likely to terminate the pro- 
cess even earlier. 

In our adopted viscosity prescription, the first stars be- 
gin to form from the same instabilities that drive angular- 
momentum loss. Of course, if the disc completely fragments 
into stars then the angular momentum transfer stops, but 
the timescale for this is likely longer than the lifetimes of 
the very first massive stars. In addition, the presence of ad- 
ditional feedback on star formation, for example fro m mag- 
netic fields, may play a role in setting this timescale llGnedin| 
l200lfl . Therefore, we assume that a more realistic estimate 
of the time available for viscosity to work is comparable 
to the timescale for collapse and evolution of the first self- 
gravitating bodies in the disc. These "first stars" should have 
very low metallicities because the haloes under consideration 
at z ~ 20 had very few progenitors massive enough to have 
formed stars before z ~ 20. Thus, the stellar evolution in 
these stars is quite rapid, with lifetimes between ^1 — 30 
Myr ; depending on many assumptions dShapiro fc Shibata 
|2002[ iRarafee et al.ll200ll: iHeger Wooslevl 120021: ISf^eTeT 
2002). We explore the range t — 1 — 30 Myr for the time 
available for viscosity to transfer angular momentum away 
from the inner disc. 

If angular momentum transfer is allowed to oper- 
ate for a timescale t, then a large fraction of mass 
within fvi s would lose its angular momentum and form 
a central, massive object. This object would be pressure- 
supported and short-lived 5 , and inevitably collapse to a 

5 It is perhaps a concern that the presence of a pressure- 
supported object that grows over the viscous timescale might in- 
fluence the temperature of the disc. We can estimate the heating 
of the disc if we assume that this objects radiates at Eddington 
luminosity with an efficiency of ~ 0.1 and that the disc is heated 
due to Thomson scattering. We find that it will take a time of 
order ~ 5t v is to increase the temperature of the disc by 50% (this 
timescale is reduced to ~ 3t v i s if the object radiates at the Ed- 
dington limit). Since these timescales are greater than the viscous 
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black hol e owing to the post-Newto nian gravitational in- 
stability dShapiro fc Teukolskvl Il983h . The exact mass of 
the resulting black hole depends on the mass loss expe- 
rienced during the short lifetime of this "star". This is 
an unknown quantity, but the mass loss is typically as- 
sumed to be i n the range 10% — 90% of the initial mass 
of the object jBarafee et al.ll200ll : iHeeer fc Wooslevl 120021: 
IShaoiro fc ShibataTT2o"o2li . Hereafter we assume that the 
mass of the resulting black hole is some fraction n of the 
total mass with zero angular momentum, and take k = 0.5 
as our fiducial value where necessary. 



4 MASS FUNCTION OF SEED BLACK HOLES 

We define the mass of the black hole Mbh as the fraction k 
of the mass within f v i s (i), i.e., the mass that lost its angular 
momentum over time t. Using equations 1131 and J7J we 
obtain for the generalized profile 

2 V 



5 x 10" 4 M 



- (JL) (AL\ 
v V0.5/ \Myi J 



3,, 



(14) 



A 



3», 



[20] 3 "- 3/2 [1.1 x 10" 



1 27J-1 



where rj = £/(3 — £). For the B01 case with £ = 1 the result 
reduces to 

«. . -,.xfe)(4r(«)";) 



_M^ 3/2 f 

0.04 



V 10 7 M e 



I Myr 



Thus, for typical values, the black hole masses are quite 
large, ~ 1O 4 M0 , and of order ~ 10~ 3 the mass of the host 
halo. Notice however that the black hole mass is quite sen- 
sitive to the value of the profile slope parameter £ (e.g. the 
last term in ea. 1141 . For example, with £ = 1.4, the black- 
hole mass is ~ 10" 4 of the halo mass. A similar calculation 
using an exponential-like profile (with £ = 2) gives negligible 
seed black hole masses, unless the viscous timescale is much 
larger than the Hubble time at that high redshift. 

In general, for a halo of a given mass M v forming a 
black hole at redshift z we expect a spread in the parame- 
ters characterizing the distribution of angular momentum. 
We compute the conditional probability distribution of seed 
black hole masses given M v at z for the £ = 1 case, where the 
distributions of the parameters A and fi are well-known. Note 
first that the black- hole mass scales as Mbh oc [fib (/x)/A] 3 , 
where we define for convenience /3 = fib(fi) and 7 = /3/A, 
such that Mbh oc 7 3//2 . The N-body simulations used by 
B01 show that the probability distributions of A and j3 are 
well-parameterized by log-normal distributions of the form, 



P(x)dx = 



1 



27ra x x 



exp 



ln 2 (x/a;o) 



2al 



dx 



(16) 



timescales we consider here (set by the timescales for the forma- 
tion and evolution of the first self-gravitating bodies in the disc) 
we conclude that the presence of this central short-lived object 
does not affect the thermal evolution of the disc. 
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Figure 5. Probability distributions of seed black hole masses 
for the threshold halo mass at 2 = 17. The assumed angular- 
momentum distribution is £ = 1, and the assumed viscosity 
timescale is 5 Myr. 



where x§ = 0.035, 1.25 and a x — 0.5, 0.3 for x = X,/3 respec- 
tively. Given these distributions, it is easy to compute -P(7), 
and from this the probability distribution of the seed black 
holes given for a halo of mass M v at redshift z for a given 
viscous timescale t. The result is 



P(Mbh|M v )dM b h 



1 



exp 



27TffM bh M bh 

ln 2 (Mbh/M b °h) 



2 < 



dMbh. (17) 



where M° h is 



M, 



3.8 x 10 4 M Q (—)(— 
VO.5/ VO.03 



3/2 



-Mv 



10 7 M 







18 



3/2 



t 



10 Myr 



(18) 



and <T Mbh =0.9. 

Figure 5 shows the seed black hole probability distri- 
bution normalized to a peak value of one, for a host halo 
of mass Af v = 5 x 1O 7 M at z = 17 and a time available 
for angular momentum transfer t = 5 Myr. In this case, 
t dyn (z = 17) w 73 Myr, tq(M = 5 x 1O 5 M ,z = 17) * 35 
Myr and therefore / ~ 0.014 yielding a mean seed black 
hole mass of Afbh ~ 3 x 1O 4 M0 with possible masses in the 
range 10 3 — 10 6 Mq . Such a distribution function of seed 
black hole masses should provide a useful starting point 
for studies aimed at modelin g the quasar properties within 
a co smologic al context IIBr omlev. Somervil le fc Fabian l 

[ 20031: iMenou. Haiman fc Narayanan! 1200 ll : 

Kauffmann fc HaehnehJ l2000t) and can be used in of the 
merger history of supermassive black ho les iMadau fc Reesl 
l200ltlVolonteri. Haardt fc Mada ull2"003h . 

Consider now black holes forming only in haloes of crit- 
ical mass MJ rlt . If we approximate M" 1 * with equation 1101 . 
then the seed black hole given in equation 114H is indepen- 
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dent of redshift of formation, host halo mass, and spin pa- 
rameter, and has a "universal" value of 



3/2 



lOMyr 



.(19) 



However, it is expected that since M^ rlt depends on the 
spin parameter A and fib (/it) there will be a spread in the 
critical host halo masses. The critical mass My scales 
as M" lt oc 7 _3,/2 . If we use the probability distribution 
for P{"i) calculated as explained above, it is straightfor- 
ward to calculate the distribution in the host halo masses. 
The result is a log-normal distribution with a mean of 
M crit,o ^ 5 x 10 7 Mq T^o [18/(1 + z)] 3/2 and a standard 
deviation of cr M crit = 0.9. Thus, if we focus on black holes 
forming within critical mass haloes only, our model predicts 
a near universal mass of ~ 1O 5 M0 , with a log-normal dis- 
tribution of the masses hosting these seeds. 



5 TOTAL MASS IN SEED BLACK HOLES 

In this section we calculate the comoving mass density of 
mass of seed black holes and make a comparison with the 
observed density of mass in black holes today in the local 
Universe. In this calculation we do not include any addi- 
tional black hole growth via luminous accretion (which must 
occur at some level in order to match the observed quasar 
population). We therefore expect the calculated density to 
be smaller than the observed value. If, for example, most 
black holes experience a luminous accretion phase charac- 
terized by a few e-folding times, we want our cumulative 
mass density to fall below the locally observed black hole 
density by roughly an order of magnitude. 

As halo masses grow with time, some haloes become 
more massive than the minimum mass threshold MJ 11 * and 
become the sites of potential seed black hole formation. Since 
the threshold mass decreases only slowly with redshift (see 
figure lO . the mass density due to seed black holes is dom- 
inated by the seeds formed at the particular redshift of in- 
terest. Naturally, all haloes more massive than the critical 
host halo mass scale were able to form a seed black hole at 
some higher redshift. However, these mass scales correspond 
to higher-sigma fluctuations and therefore do not contribute 
significantly to the mass density. In order to include only the 
contribution coming from seed black holes formed in haloes 
that have just crossed the minimum mass threshold, we com- 
pare the typical growth rate of halo mass with the growth 
rate of MJ rit . We find that the probability of a halo crossing 
the threshold having previously been massive enough to host 
black hole forming progenitors is vanishingly small (~ 10 -6 ). 
Therefore it is safe to assume that the contribution to the 
seed black hole mass density at each redshift comes predom- 
inantly from those seeds formed in haloes which have just 
crossed the minimum mass threshold near that particular 
redshift. 

The comoving density of mass in black holes as a 
function of redshift can be determined by integrating over 
the formation rate of haloes at the threshold mass from 
some initial redshift z; to the redshift of interest z, i.e., 
ph\i{z) — J (dp/dz)dz. We approximate the evolution in 
the co-moving number density of seed black holes, as 
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Figure 6. Total mass density in seed black holes as a function of 
redshift. Thick lines represent our fiducial £ = 1 (B01) angular- 
momentum distribution for various viscous timescales, while the 
thin line corresponds to § = 1.3 and a viscous timescale of 20 
Myr. The shaded area represents the observed range of total mass 
density in black holes at z = 0. The WMAP reionization range is 
marked by the vertical dotted lines. 



dp{z) 
dz 



d_ 

dz 



M$ rlt (z) 



dn(z) 
dM 



M hh (M,z)dM 



(20) 



Here [dn(z)/dM]dM is the number density of haloes with 
mass in the interval between M and M + dM at redshift z, 
given by 



^ dM = £l fM M), z] 
dM M 



da 2 (M) 



dM 



dM 



(21) 



and <r(M) is the mean square fluctuation amplitude at a 
mass scale M. Here /i is the fraction of mass associated 
with haloes of mass M corresponding to the given range 
in a(M) at redshift z. We use the Press-Schechter (1974) 
approximation for f\ (see, e.g. lLacev fc Coldll99&t) . Since 
the value of M^ Ilt (z) is associated with increasingly higher 
sigma fluctuations at high redshifts, our results are not sen- 
sitive to the choice of z%. We choose z; = 30, but our results 
are unchanged for Zi — 50. 

Figure 6 shows the comoving density of seed black holes 
obtained by integrating eq. 1201 . As expected, this function 
is strongly dependent on the time available for viscosity to 
redistribute angular momentum. At high z there is a rapid 
growth, which flattens off toward lower redshifts. The hor- 
izontal band marks the ob servationally - inferred density of 
mass i n black holes today |s7) ltan||l982l: jCliolcslii_j£_Tjxrnei 
1992[ iKormendv fc Richstond ll99St Fabian & Iwasawa 



199a ISalucci et alJ Il999t iMerritt fc Ferraresd |20(L 
^i^r^IVelmdndn2002l) . The total mass in black holes 
increases with time, until the cool gas supply stops due to 
some feedback mechanism that prevents and/or destroys 
the production of molecular hydrogen. Here we assume 
that cosmological reionization at z TC can be responsible for 
preventing the further cooling and collapse of haloes of 
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mass M£ for redshifts z < z Ic . This radiative feedback is 
expected to be especially effective in the relatively small 
systems under consideration. The ±1<t range for the red- 
shift of reionization as measured by WMAP (z IC = 17 ± 5, 
ISpereel et all 120031) is marked by the vertical lines. We 
expect the model prediction for the seed mass density at 
z la to be roughly an order of magnitude below the accepted 
observed value today, in order to allow for the additional 
growth by accretion during an AGN phase. We see that 
for z IC = 17, say, this requirement is fulfilled once viscosity 
operates for a timescale of 1 ^ t ^ 5Myr. 

Also shown in Fig. 6 is the case £ = 1.3 with a viscous 
timescale of t = 20 Myr. The curve is flatter in this case be- 
cause M^ rlt for this profile grows faster than the £ = 1 case 
(see figure 4 ) , and as a result the contribution from new seed 
black holes is smaller. In addition, the mass of the black hole 
associated with each halo above M" 1 * scales as oc (1 + z) 3,8 
for £ = 1.3 compared with oc (1 + z) 1 ' 5 for £ = 1 (see equa- 
tion ll5> . We see that the model predictions for £ = 1.3 are 
acceptable only for long star-formation times, t ^ 20 Myr. 
Flatter profiles, with £ significantly larger than 1.3, either 
require a very long viscous timescale for seed formation, or 
they fail to produce enough mass in black holes by the mech- 
anism proposed here. In this case one has only small early 
seeds and therefore must accrete for many e-folding times. 



6 BLACK HOLES AND SPHEROIDS 

Our model provides a prediction for massive seed black holes 
with a characteristic linear correlation between the masses 
of the seed black holes and their host haloes. In this section 
we show how our model for high-redshift seeds may lead nat- 
urally to a relation between the black hole mass and stellar 
spheroid mass at z = by following their evolution through 
halo merger trees. We associate a galaxy's disc mass with 
gas that is accreted as diffuse material since the last major 
merger event. Each incoming halo is assumed to host a mass 
of cold gas in its center, which can contribute to a spheroid 
component if the merger occurs on a short dynamical fric- 
tion timescale. Similarly, if the merged halo hosts a black 
hole seed (see below) and reaches the central object in time, 
then this seed can merge with the central black hole. In this 
fiducial example, we assume that when black holes merge, 
their masses add linearly. 

We gen erate halo merger tre e s usin g the algorithm de- 
veloped by ISomer ville fc Kolattl <1999h based on the ex- 
tended Press-Schechter formalism l|Press fc Schechteilll973 
iBond et alJll99ll:lLacev fc Colelll993) . According to this al- 
gorithm, the probability for a mass increase by AM in a 
time step At is 

P{Aa 2 , AS c )d{Aa 2 ) = ^ c ■ e ~ M l/^" 2 d(Aa 2 ), (22) 
\>2nAa 3 

where A8 C = 5 c (t)—8 c (t+At) is the linear fluctuation growth 
according to the chosen cosmological model, corresponding 
to the time step At, and Acr 2 = a 2 (M)-a 2 (M + AM) is the 
increment in mean-square fluctuation amplitude according 
to the given linear power spectrum, corresponding to the 
mass increment AA1. We choose the time step such that 

AS C < Jm~^^- (23) 



for a chosen minimum halo mass M m i n ; a mass increase by 
AM < Mmi n is considered to be a smooth accretion process 
rather than a merger event. The merger tree is generated 
by randomly selecting progenitor masses according to the 
above probability distribution at each discrete time. The 
parent halo is identified with the most massive progenitor 
at each time while the rest of the mass is either in subhaloes 
or diffuse mass. The tree is truncated at a high redshift 
when the mass of the major progenitor is less than M m i n . 
We run 50 Monte Carlo realizations of the merging histories 
for today's haloes in the mass range 10 11 — 1O 14 M0 , with 
M min = 1O 6 M . 

For each merger, we take into account the finite time for 
the incoming subhalo to spiral in due to dynamical friction 
and reach the center of the host halo. We model each halo 
as a singular isothermal sphere with a radius and velocity 
given by the virial relations discussed i n §2. The dynamical- 
friction time can then b e written as feinnev fc Tremaind 
ll987tlLacev fc Coldll993T) . 



tdf = 0.42 Gyr 




where R is defined as the ratio of the mass of the major pro- 
genitor to the mass of the incoming halo. For our calculation 
we approximate the Coulomb logarithm as InA ~ \nR. The 
function g(e,ri) = e°' 5 r/ 2 characterizes the e ffect that non- 
circular orbits have on the orbital decay (e.g. lTormenlll99H : 
iGhigna et a!] 119981. Here e characterizes the circularity of 
the orbit, and is defined as the ratio of the angular momen- 
tum of the orbit to the angular momentum of a circular orbit 
with the same energy. The parameter rj characterizes the en- 
ergy of the orbit by defining te = f)R v to be the radius of a 
circular orbit with the same energy as the actual orbit. We 
will adopt the typical parame ters r\ = 0.6 and e = 0.5 (see, 
e.g.. lZentner fc Bulloc3l2003ft . 

Material that ends up in a disc is assumed to be the 
gas that comes in as diffuse accretion since the last ma- 
jor merger. The last major merger is defined as a merger 
with an incoming object that is at least 50% the mass of 
the host. The diffuse gas may partly be the gas lost from 
dark haloes smaller than M~2x 10 s [(1 + z)/18]~ 3/2 due 
to p hoto-evaporation starting at the epoch of reionization 
(e.g. iBullock et alJ l200cl : IShaviv fc Dekell 120041) . Although 
the gas escapes from the small haloes, it may remain bound 
to the larger host halo, and we assume that it is accreted 
onto the disc at about the same time the low-mass haloes 
are accreted. For haloes less massive than ~ 3 x 10 n M o , 
which roughly coincides with the upper limit for disc for- 
mation at high redshift, we indeed expect short infall times 
because of short cooling times and negligible shock heating 
jBirnboim fc Dekellliooa) . 

The spheroid mass is assumed to be the straightforward 
sum of the cold baryons in the merged subhaloes throughout 
the merger history that sank to the center of the subhalo. 
We take the spheroidal component to be in place by the 
redshift of the last major merger assuming that a major 
merger would disrupt the system into an irregular galaxy. If 
the estimated sinking time of the incoming subhalo is shorter 
than the time interval between the subhalo merger and the 
last major merger we add to the bulge a mass foM s , where 
fo — 0.03 is the assumed cold-gas fraction, and M s is the 
mass of the merged halo. The choice of the mass fraction 
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in the definition of the major merger makes only a small 
difference to the final slope of the relation between black 
hole and spheroid mass, because the merger history tends 
to converge to the final mass after the last major merger. 

Among the subhaloes that reach the center before the 
redshift of the last major merger, we determine if they have 
a black hole, and if so, we assume it merges with the central 
black hole over a short timescale. In principle, we could com- 
pute the black hole mass contained in each merged halo by 
following its full merger tree back in time, until some early 
redshift z> z rc . However, since our aim is simply to investi- 
gate the ramifications of our scenario, we have chosen only 
to follow the merger tree of the most massive progenitor, 
and estimate the formation times of each merged subunit in 
a statistical fashion. This choice saves a significant amount 
of computational time. 

Specifically, since seed black holes stopped forming at 
z = z rc , we can estimate the total black hole mass present 
in a subhalo of mass Ms that merges into the main host at a 
redshift z meI by the following approximation. If z mer < z Ie , 
then the black-hole mass populating the merging subhalo 
is determined by integrating over the mass function of its 
progenitors at z rc , 

f°° M 
M hh = / ^M bh {M,z TC ) 

J M^Hzro ) M 

x f 2 [a(M),z rc \a(M B ), 
da 2 (M) 



dM 



\dM. 



(25) 



Here, M is the mass of progenitors. The function 
/2[<r(M), Zro | ff(M s ), Zmer] represents the fraction of mass in 
haloes of mass M s at z mcr that wa s in haloes of mass M 
at redshift z rc (see equation 2.15 of lLacev fc Colel lll993j) ). 
The mass Mbh(M, z ro ) is the mass of the seed black hole 
residing in a halo of mass M at z ro according to eq. 11511 . If 
the merger occurs before z rc , then the black-hole residing in 
the merging subhalo is assumed to have a mass according to 
eq. 1151 . given the halo mass M B at z mcr . We only need to 
follow the merger history back to the time when the major 
progenitor is as small as M^ rlt (z), because when a progenitor 
is smaller than M" 1 *, no seed black holes will form. 

The integral above gives only an average black hole mass 
for an ensemble of host haloes of mass M s . So the approx- 
imation is valid if Mbh 3> Mo, where Mo is the seed mass 
associated with critical-mass haloes at z ro . However, for low- 
mass subhaloes, the approximation can break down because 
it allows solutions with black hole masses smaller than the 
initial seed masses, Mbh < Mo (~6x 1O 4 M0 , see Eq. 18). 
This is unphysical, because we assume that each black hole 
is built up by a series of mergers with seeds of mass Mo- We 
attempt to remedy this shortcoming in cases by randomly 
re-assigning the value of Mbh in these cases. Therefore if 
Mbh < Mo we choose a random uniform deviate r, and set 
Mbh = if r > Mbh/Mo and M b h = Mo otherwise. 

Figure 7 shows the resulting black-hole mass as a func- 
tion of the host spheroid mass from the merger tree realiza- 
tions described above (open symbols), in comparison with 
the observed data (solid points with errors, see caption). We 
have adopted model parameters of t = 10 Myr and z rc = 15. 
The model predictions in the top panel refers to a given fixed 
halo mass today, M = 3.2 x 1O 12 M0 , where only the ran- 
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Figure 7. Today's black-hole mass ver sus host-spheroid mass . 
The filled symbols are based on data from lTremaine et alJ<2002ft . 
where the bulge mass is derived from the B absolute magni- 
tude assuming the constant mass-to-light ratio advocated by 
Ivan der Marell jl99ll) . The model predictions are based on merger 
growth only, ignoring accretion. Top: for a fixed halo mass to- 
day, M v = 3.2 X 10 12 Mq , using 50 random realizations of the 
merger tree. The variations in merger history by themselves tend 
to spread the masses along the observed relation. Bottom: for 
a range of present halo masses, M v = 10 11 — 10 13 ' 5 Mq , with 
logarithmic half order of magnitude steps. The slope of he pre- 
dicted relation is in agreement with the observed slope. In cases 
where the black-hole mass is calculated to be zero, we have placed 
the symbol at the bottom of the panel for clarity. The predicted 
black-hole masses could potentially increase by about an order of 
magnitude due to accretion in an AGN phase. 



dom merger history is varied from point to point. Though 
not shown, the associated disc masses for these haloes range 
from 10 10 — 10 n MQ with no clear correlation between the 
disk and black hole masses. 

We find that the black hole mass correlates with the 
mass of the spheroid component, as observed, even for a 
fixed-mass host halo. This result arises because subhaloes 
that are massive enough to sink to the center of the parent 
halo and contribute to a bulge, are generally massive enough 
to host black holes as well. Haloes which experienced many 
mergers with small subhaloes tend to have a small spheroids 
and small black holes. This would naturally explain why 
black hole mass is observed to correlate w ith spheroid mass 
and not disc mass llGebhardt et alJl2001f) . since disc mass 
does not correlate well with the merger history in our pic- 
ture. 

The bottom panel of Fig. 7 shows several ensemble re- 
alizations for a series of host halo masses with z = values 



of 10 1 



10 Mq . The computed relation follows closely 



Mbh oc M sp , and is qualitatively similar to the observed re- 
lation at z = 0. The normalization (set primarily in this 
case by setting t = lOMyr) allows for an additional growth 
by a few e-folds during an AGN accretion phase, as desired. 
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In addition, for the lowest-mass spheroids, we see in some 
cases the absence of any black hole mass (points along the 
bottom edge represent Mbh = 0). This comes about because 
low mass host haloes are less likely to have had progeni- 
tors more massive than M" 1 * at z = z Ic . Black holes cannot 
be smaller than the characteristic seed mass ~ 6 x 1O 4 M 
while there is no fixed mass threshold for a progenitor to 
contain gas that will eventually end up in a spheroid. This 
gives rise to a slightly steeper slope in the bottom left of the 
relation. The apparent increase in the number of low mass 
spheroids with no black holes also suggests the presence of a 
lower mass limit on spheroids hosting black holes, of roughly 
M sp ~ 10 8 M Q . This is a potentially verifiable prediction of 
the model. 

It is of interest to investigate the robustness of the de- 
rived relation between Mbh and M ap to specific assumptions 
of the model presented here. For example, the result is prac- 
tically independent of the finite duration of a merger due 
to dynamical friction. Major mergers are quick and have 
an immediate effect on the growth rate of the spheroid as 
well as the black hole, while minor mergers take a long time 
to complete and therefore have a delayed effect that may 
not materialize in a Hubble time JZentner fc Bullocklliooa 
iTavlor fc Babulll200ll l2003ft . When we repeat the calcula- 
tion under the assumption that all mergers are instanta- 
neous, ignoring the finite duration of the spiral-in process, 
we find that the predicted relation between black holes and 
spheroids remains almost unchanged, with only a slight shift 
of typically less than 10% towards larger spheroid masses for 
the same black hole mass. This small shift is to be expected 
since the less massive subhaloes that now contribute to the 
spheroid mass tend to be below the critical mass for seed 
formation and thus do not add to the black hole mass. 

We find the predicted slope of the relationship to be 
insensitive to the actual value of M£ (as long as M v ^> 
M™ ). The normalization of the relation, namely the black- 
hole mass for a given spheroid mass, is roughly inversely 
proportional to M" 1 '. This is because for a lower value of 
MJ rl more haloes are able to form a seed and contribute to 
the merged black-hole mass without affecting the spheroid 
mass. 

We also find that the slope is insensitive to the value 
of z IC (as long as z rc ^ 6), while the normalization does de- 
pend on it. For a lower z le , seed formation takes place in a 
larger fraction of the haloes that eventually merge to form 
spheroids. By lowering z la from 15 to 6 we obtain an or- 
der of magnitude increase in black-hole mass. A comparable 
decrease in black-hole mass is obtained when increasing z TO 
from 15 to 20. 

However, the slope of the relation is found to be mildly 
sensitive to the original scaling between the seed black hole 
mass and its host halo mass. If we replace the predicted cor- 
relation of our model (Mbh oc M v ) with an ad hoc seed 
assignment to critical- mass haloes at high redshift, e.g., 
M b h oc (M" jt ) 3 , we find that the linear M b h oc M sp relation 

3/2 

steepens to roughly Mbh oc M sp . This is because the ra- 
tio Afbh/Afv is now monotonically increasing with M v . This 
mild sensitivity means that the predicted correlation of our 
physical model for seeds is important, but it does not have 
to be exactly Mbh oc M v . 

Finally, we point out that our result is quite sensitive 



to the way we assign mass M to a black hole which results 
from the merger of two black holes, of masses Mi and M%. In 
our fiducial calculation we have assumed that the black hole 
masses sum up linearly, M = Mi + M2. On the other hand, 
if entro py is conse r ved, the mass o f the merger product is 
M = yjMl + Mf jHawkindll971allrl . This should serve as 
a lower bound to the possible mass of the merger product, 
and should yield the flattest-possible Mbh — M sp relation. 
Consider a black hole in the center of the halo formed from 
iV m mergers of seed black holes of mass Mbh ~ 10 5 Mq . If 
entropy is conserved, and -/V m ^> 1 then the black hole mass 
would be Mbh = Af sccd (l + A m ) 1/ ' 2 oc Nl[ 2 . In our model, 
spheroids grow via a series of mergers with the gas from 
nearly equal-mass haloes, roughly as Af sp ~ /oM s iV m oc iV m , 

1 /2 

so we expect Mbh oc M sp in this limit. When we follow this 
recipe for mass adding using our full merger-tree method, 
we indeed obtain a result very close to our expectation: 
Mbh oc Af s ° p 55 . Thus, in this limit, the derived relation would 
be somewhat flatter than observed. In order to comply with 
the observed relation, one would have to appeal to a sce- 
nario where black-hole growth via accretion is more eff icient 
in hi gh-mass spheroids (e.g., IVolonteri. Haardt fc Madaul 
120031) . 

With the above tests in mind, we conclude that the de- 
rived linear relation between black hole and spheroid mass 
depends primarily on the fact that the seeds were massive, 
that they formed at high redshift, and that their masses add 
roughly linearly when they merge. In this case the build-up 
of black holes and spheroids is occurring in concert through 
mergers, thus setting the linear relationship. A similar linear 
relation would hold for the final black-hole masses provided 
that it has not been changed by the subsequent growth ex- 
perienced by black holes via accretion. If the accretion is 
proportional to the mass of the seed, as would be the case 
for accretion at the Eddington rate, then we expect the slope 
of the relation to persist. Of course, the effect of accretion 
depends on what fraction of the black hole mass was ac- 
creted. If the increase in mass is only by a factor of a few, 
then one could assume that the role of mergers in setting 
the observed relationship is important. On the other hand, 
if the growth by accretion is by more than an order of mag- 
nitude, then the accretion process would have the dominant 
role in determining the slope of the relation. According to 
our model, the normalization of the observed relationship is 
determined additively by two free parameters: the viscous 
timescale and the growth via accretion. When adopting a 
physically motivated value for the viscous timescale, we find 
that growth via accretion is limited to be by an order of mag- 
nitude or less, thus suggesting that mergers do play an im- 
portant role in setting the relationship. We note that in mod- 
els where the seeds are of smaller masses, ~ 100 Mq (e.g., 
IVolonteri. Haardt fc Madaul 12003ft . the growth by mergers 
is negligible and the accretion process has to be tuned with 
a specific feedback scenario in order to obtain the correct 
relation between black holes and spheroids. 

Although we focused here on black hole and spheroid 
growth, we do not expect all the seed black holes accreted 
by a dark matter halo to be incorporated in the spheroid. 
A number of black holes could be orbiting in haloes of 
big galaxies. They could have arisen from subhaloes that 
were disrupted after ent ering the hos t halo. Such a scenario 
has been investigated bv llslam etafl j2003al) and the possi- 
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ble obs ervational consequences are explored in llslam et al. 




7 CONCLUSION 

We presented a physical model for the production of massive 
seed black holes at high redshifts, with a characteristic mass 
~ 10 Mq . We argued that this model provides a useful 
starting point for explaining the observed relation between 
the masses of today's black holes and the spheroidal compo- 
nents of galaxies. Massive seeds may also help explain the 
existence of the supermassive black holes associated with 
luminous AGN at z> 6. 

Our model for seed formation is based on the general 
idea that the mass in every early-collapsing halo is expected 
to have a distribution of specific angular momentum, with at 
least some fraction of its gas having a much smaller specific 
angular momentum than is characterized by the global spin 
parameter of the halo. After the low-spin gas cools, it can fall 
into the central region and form a dense disc, where certain 
viscous processes become effective and may potentially lead 
to black hole formation. 

In our exploration of this idea, we assumed that the 
specific angular momentum of the gas is similar to that 
measured for the dark-matter in N-body simulations. In this 
case, if angular momentum is conserved during the gas infall, 
the central gas discs become self-gravitating and Toomre un- 
stable in all the haloes above a critical mass. This critical 
mass has a median value of ~ 7 x 1O 7 M0 at z ~ 15, and it 
is expected to be scattered log-normally with a In standard 
deviation of ~ 0.9. The sites of seed production are thus 
much more massive than the typical collapsing mass at that 
time, and are therefore associated with high-sigma peaks 
in the fluctuation distribution. Gravitational instabilities in 
the disc drive an effective viscosity that redistributes angu- 
lar momentum. As a result, material from the low angular 
momentum tail of the distribution loses its angular momen- 
tum on a viscous timescale and contribute to the formation 
of a black hole. We associate this viscous timescale with the 
timescale for disruption and/or heating of the disc due to 
substantial star formation. We expect this timescale to be 
on the order of 1 — 30 Myr, but its exact value is not well 
determined. 

We found that the typical seed black holes, which 
form in critical mass haloes, have a characteristic mass of 
~ 10 s Mq , quite insensitive to the epoch of formation. This 
predicts a lower bound to possible masses of massive black 
holes, which can be confronted with observations today. 

This scenario should operate continuously as additional 
haloes grow to above the critical mass and create new poten- 
tial sites for black-hole formation. The process is expected 
to stop, however, when the Universe becomes sufficiently 
reionized that H2 molecules are destroyed and cannot pro- 
vide efficient cooling. We therefore expect seed black holes 
to form before z re ~ 13 — 21. 

We showed that the comoving density of mass in seed 
black holes grows gradually with time, in association with 
the birthrate of critical mass haloes. This density becomes 
comparable to the observed mass density of supermassive 
black holes by z ~ 15, near the redshift of reionization. 
Depending on the adopted viscous timescale, there is room 



for an order of magnitude growth in black hole masses via 
quasar powering accretion. We find that a viscous timescale 
of t ~ 5 — 10 Myr provides quite reasonable results. 

We explored the development of a correlation between 
black-hole mass and spheroid mass in today's galaxies by 
tracking our early forming massive seed population through 
a halo merger tree. Under the assumption that spheroids 
form via major mergers, we showed that the merger process 
leads naturally to a linear relation between black-hole mass 
and spheroid mass. This result is mainly due to the forma- 
tion of massive seeds only in rare, high-peak haloes at high 
redshifts, and is not too sensitive to the details of our seed 
formation model. Starting with ~ 1 — lOOM© seeds would 
not lead to a linear correlation via mergers alone; additional 
feedback associated with AGN accretion is needed. The seed 
black holes in our model have masses that are proportional 
to the virial masses of their host haloes at formation. We 
demonstrated that this linear relation is preserved if black 
holes merge in concert with the hierarchical growth of haloes 
and spheroids. 

One potentially interesting ramification of our scenario 
is that a population of massive black holes at high redshift 
might serve to power small quasars, fueled perhaps by low- 
angular-momentum gas liberated by merger events tha t oc- 
cur after the seeds are in place. IMadau et al.l (|2003j) dis- 
cussed how a population of "mini-quasars" like this might 
be responsible for reionizing the universe at z IC ~ 15. How- 
ever, in their model, the mini-quasar power-sources were 
intermediate-mass ~ lOOM© black holes. Under similar as- 
sumptions, the massive seed black holes of our model would 
produce even more ionizing photons per site, and perhaps 
lead more naturally to an early reionization. In our cur- 
rent scenario the reionization epoch was an input param- 
eter detached from the black-hole formation process, but, 
in principle, the accretion onto the seeds themselves could 
be responsible for the reionization that eventually inhibits 
further black-hole formation. 

Our quantitative predictions should not be interpreted 
too literally. For example, our fiducial model assumes that 
the initial angular-momentum distribution of the gas mir- 
rors that of a typical CDM halo, and that it is preserved 
during disk formation. There is no evidence in local galax- 
ies for discs with inner densities as high as invoked by our 
model for proto-galaxies. This might be explained by pro- 
cesses that have altered the gas angular-momentum distri- 
bution preferentially at late epochs, such as feedback effects 
(e.g. iMaller fc Dekell2002l) . but one should clearly assign 
a substantial degree of uncertainty to our model assump- 
tion. In order to address this uncertainty, at least qualita- 
tively, we explored the effect of varying the initial angular 
momentum distribution. We found that if the low angular 
momentum tail is as minor as in an exponential disc of a 
constant-density core (£ = 2), the resulting seed black holes 
are too small to be of relevance. An intermediate case, with 
£ = 1.3, can produce sensible black-hole seeds, compara- 
ble to our fiducial case (£ = 1), provided that the viscosity 
timescale is slightly longer, 20 ^ t ^ 30 Myr. For £ = 1.3, 
the critical halo mass becomes ~ 1O 8 M0 at z = 15, and the 
inner ~ 10% of the disc mass manages to lose its angular 
momentum by viscosity and produce a black hole of mass 
~ 1O 6 M0 . This seed mass is larger than in the fiducial case, 
but the seeds in this case are rarer events (see figure 13. II . 
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Therefore, in order to recover a similar total mass density 
in black holes, the viscosity must operate for a longer du- 
ration. The actual low angular-momentum tail of the gas 
distribution remains an unknown, which will hopefully be 
constrained by future, more realistic hydrodynamical sim- 
ulations, properly including feedback and other important 
physical effects. 

In conclusion, while the model presented here is clearly 
idealized, it does provide a useful example for the kind of 
mechanism needed to connect cosmological structure for- 
mation on galactic scales to black-hole formation on much 
smaller scales. Observations suggest the presence of massive 
black holes at high redshifts and therefore the need for the 
formation of massive seed black holes at still earlier times. 
According to our scenario, the observed correlations of black 
hole and galaxy properties are directly linked to the hierar- 
chical nature of structure formation itself, and owe their ex- 
istence to rare systems forming in association with the first 
light. 
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